Scalable Parallel RK Solvers for ODEs Derived by the Method of Lines
نویسندگان
چکیده
This paper describes how the specific access structure of the Brusselator equation, a typical example for ordinary differential equations (ODEs) derived by the method of lines, can be exploited to obtain scalable distributed-memory implementations of explicit Runge-Kutta (RK) solvers. These implementations need less communication and therefore achieve better speed-ups than general explicit RK implementations. Particularly, we consider implementations based on a pipelining computation scheme leading to an improved locality behavior.
منابع مشابه
Parallel Low-Storage Runge - Kutta Solvers for ODE Systems with Limited Access Distance
We consider the solution of initial value problems (IVPs) of large systems of ordinary differential equations (ODEs) for which memory space requirements determine the choice of the integration method. In particular, we discuss the space-efficient sequential and parallel implementation of embedded Runge–Kutta (RK) methods. We focus on the exploitation of a special structure of commonly appearing...
متن کاملExploiting Limited Access Distance of Ode Systems for Parallelism and Locality in Explicit Methods∗
The solution of initial value problems of large systems of ordinary differential equations (ODEs) is computationally intensive and demands for efficient parallel solution techniques that take into account the complex architectures of modern parallel computer systems. This article discusses implementation techniques suitable for ODE systems with a special coupling structure, called limited acces...
متن کاملParallel Iterated Runge Kutta Methods and Applications
The iterated Runge Kutta IRK method is an iteration scheme for the numerical solu tion of initial value problems IVP of ordinary di erential equations ODEs that is based on a predictor corrector method with an Runge Kutta RK method as corrector Embed ded approximation formulae are used to control the stepsize We present di erent parallel algorithms of the IRK method on distributed memory multip...
متن کاملDefect Sampling in Global Error Estimation for ODEs and Method-Of-Lines PDEs Using Adjoint Methods
The importance of good estimates of the defect in the numerical solution of initial value problem ordinary differential equations is considered in the context of global error estimation by using adjoint-equation based methods. In the case of solvers based on the fixed leading coefficient backward differentiation formulae, the quality of defect estimates is shown to play a major role in the reli...
متن کاملOn second derivative 3-stage Hermite--Birkhoff--Obrechkoff methods for stiff ODEs: A-stable up to order 10 with variable stepsize
Variable-step (VS) second derivative $k$-step $3$-stage Hermite--Birkhoff--Obrechkoff (HBO) methods of order $p=(k+3)$, denoted by HBO$(p)$ are constructed as a combination of linear $k$-step methods of order $(p-2)$ and a second derivative two-step diagonally implicit $3$-stage Hermite--Birkhoff method of order 5 (DIHB5) for solving stiff ordinary differential equations. The main reason for co...
متن کامل